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DEVELOPMENT OF DATA PROCESSING, INTERPRETATION AND 
ANALYSIS SYSTEM FOR THE REMOTE SENSING OF 
TRACE ATMOSPHERIC GAS SPECIES 

By 

Joseph C. Casas, 1 Mary S. Saylor, 2 and Earl C. Kindle^ 

INTRODUCTION 

The Cooperative Agreement, NCCl-34, between the NASA/Langley Research 
Center (LaRC) and Old Dominion University (ODU), entitled "Development of 
Data Processing, Interpretation and Analysis System for Remote Sensing of 
Trace Atmospheric Gas Species," represented a major research effort since 
its award in January 1980. Research done under this Cooperative Agreement 
was a followup to work performed under NASA grants NSG-1127 and NSG-1395, 
"Data Reduction Analysis and Application Technique Development for Atmo- 
spheric Trace Gas Constituents Derived from Remote Sensors on Satellite or 
Airborne Platforms." A progress report describing initial accomplishments 
under this Cooperative Agreement was submitted to the Technical Monitor in 
June 1982. 

The major emphasis of this investigative research was the advancement 
of remote sensing technology. In particular, the gas filter correlation 
radiometer (GFCR) technique was applied to the measurement of trace gas 
species, such as carbon monoxide (CO), from airborne and Earth orbiting 
platforms. Through a series of low altitude aircraft flights (1 km to 6 
km), high altitude aircraft flights (8 km to 12 km), and orbiting space 
platform flights, data were collected and analyzed, culminating in the first 

^Sperry Corporation, Technology Applications Inc., NASA/Langley Research 
Center, Hampton, VA 23665. 

^Research Associate, Department of Geological Sciences, Old Dominion Uni- 
versity, Norfolk, VA 23508. 

^Eminent Professor, Department of Geological Sciences, Old Dominion Uni- 
versity, Norfolk, VA 23508. 


global map of carbon monoxide concentration in the middle troposphere and 
lower stratosphere. 

Research personnel had significant responsibility in four major areas 
of this remote sensing program known as the Measurement of Air Pollution 
from Satellites (MAPS) experiment. These areas were: 

(1) data acquisition 

(2) data processing, analysis, and interpretation algorithms 

(3) data display techniques, and 

(4) information processing. 

This report contains an overview of the contributions of the researchers in 
each of these areas, under this research project. 

DATA ACQUISITION 

As described in the Progress Report of May 1982, a prototype of the 
space-borne MAPS instrument was flown in aircraft flight tests. A Cessna 
402B carried the MAPS brassboard GFCR in low altitude tests conducted over 
the St. Petersburg/Clearwater, Florida area in August 1978. As part of the 
MONsoon Experiment (MONEX) in the summer of 1979, mid-altitude measurements 
were made from a Convair 990, primarily in the vicinity of the Indian Ocean. 
Data acquired in these programs, in addition to satisfying the primary ob- 
jectives of obtaining localized air pollution measurements, formed a basis 
for comparison with subsequent measurements of the Earth orbiting MAPS 
sensor. 

In November 1981, the space shuttle Columbia carried into earth orbit 
its first scientific payload. MAPS was one of the five experiments com- 
prising this payload sponsored by NASA's Office of Space and Terrestrial 
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Applications (OSTA). The MAPS experiment acquired thirty-five hours of 
nadir viewing data and underwent five in-flight calibrations. The OSTA 
payload was launched again in October 1984 on board the Space Transportation 
System (STS) orbiter Challenger. During this 8-day mission, the MAPS flight 
tape was filled to capacity with 105 hours of measurements, selectively 
recorded according to spacecraft attitude and underlying Earth targets. 
Figure 1 shows the data-taking ground tracks of the 41-610STA-3 flight. 

ODU researchers were an integral part of the team which supported these 
two shuttle flights of the MAPS experiment. They participated in pre-flight 
hardware installation and functional testing, in mission simulations, and in 
pre-launch and launch support, all at the Kennedy Space Center (KSC) . They 
traveled to the Johnson Space Center (JSC) to participate in both simulated 
and actual mission activity of the Payload Operations Control Center (POCC) . 
Landing support, including recovery of flight magnetic tape and film and de- 
integration of hardware, was supported at the Dryden Researcn Center (DRC) 
for the OSTA-1 mission and at KSC for OSTA-3. 

One aspect of the MAPS experiment data validation program was a series 
of aircraft underflights of selected STS orbits. Correlative measurements 
were obtained with a duplicate sensor and in situ instrunentation. For the 
OSTA-1 mission, aircraft operated over four locations: Farmville, Virginia; 

Key West, Florida; the Pacific Ocean, 600 miles west of San Francisco; and 
Melborne, Australia. ODU personnel assisted in the installation and check- 
out of the experiment hardware on the NASA Lear jet which overflew the two 
eastern United States test sites. Working as experiment crew, ODU 
researchers operated the hardware and gathered meteorological data and air 
samples for gas chromatographic analysis. 
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d tracks of the 


For the correlative measurement programs of both shuttle missions, 
meteorological analysis and weather forecasting services were required to 
determine the prime locations for the aircraft underflights. A direct tap 
to the Satellite Field Services Office in Camp Springs, Maryland, was set up 
in the University's weather center during the OSTA-1 mission. Images from 
weather satellites were transmitted real-time via this data link. These 
satellite facsimiles, along with data routinely collected and catalogued by 
the weather center, enabled the selection of those underflight targets 
within the fuel range of the aircraft with greatest likelihood of essen- 
tially cloud-free skies. 

In addition to meteorological paraneters, available real-time, upper 
air analyses of temperature, dewpoint depression, and wind field data were 
required to construct the atmospheric models needed to reduce the global 
data sets. These analyses were ordered from the Navy's Fleet Nunerical 
Oceanography Center (FNOC) in Monterey, California, and from NOAA's National 
Climatic Center (NCC) in Asheville, North Carolina. The computer compatible 
tapes transmitted were decoded and converted to locally accessible computer 
files. The polar stereographic and rectangular gridded analyses were vali- 
dated through intercomparisons and through comparisons to real-time meteo- 
rological measurements and climatology. 

Other support data were required. JSC supplied the Best Estimated 
Trajectory (BET) product for the OSTA-1 mission and the on-orbit Postflight 
Attitude and Trajectory History (PATH) data for the OSTA-3 mission. World- 
wide radiosonde launch information was obtained from the NMC. Ohio State 
University supplied a grid 1* x 1“ of mean topographic elevations or bathy- 
metric depths. The primary source of this latter data set was the Defense 
Mapping Agency Aerospace Center in St. Louis. As was true with the upper 
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air analyses, these data sets were also received on computer compatible 
tapes, each having a unique storage format which required conversion by ODU 
personnel for access and use at LaRC. 

DATA PROCESSING, ANALYSIS, AND INTERPRETATION ALGORITHMS 

Extensive testing of the MAPS instrument was conducted prior to the 
OSTA-1 launch. A standardized acceptance test procedure was followed before 
and after every relocation of the instrument to evaluate functional perform- 
ance. To verify that the MAPS hardware would function properly in the 
environment of space flight, the instrument underwent a series of compliance 
tests; vibrational, thermal vacuum, and electromagnetic interface (EMI). 
After installation of the scientific hardware on the OSTA pallet, integra- 
tion testing and a final end-to-end check out were conducted in the Vehicle 
Assembly Building at KSC. 


Instrument Calibration 

Prior to the STS-2/0STA-1 flight, the instrument was calibrated in a 
laboratory at room temperature. Following the flight, it was calibrated in 
a thermal vacuum chamber at a variety of instrument temperatures. During 
the actual mission, data were acquired while the instrument was viewing deep 
space which allowed the determination of the zero position. 

Preflight 

The MAPS instrument preflight calibration tests were carried out in a 
laboratory at an ambient temperature of 298 K. These tests determined the 
magnitude of the output noise, the radiometric responsivity, the degree of 
linearity of the output signal response, voltage offsets, and the repeat- 
ability of the instrument for each of its three channels: V, containing an 
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evacuated cell; AV, containing a high pressure (266 torr) CO cell; and AV , 
containing a low pressure (76 torr) CO cell. 

The specially designed calibration system controls independently the 
composition, total pressure, and temperature of the gas mixture contained in 
a calibration cell of length 0.5 meter. The system also varies the tempera- 
ture of a high quality calibration blackbody between 238 K and 345 K. The 
measured temperatures and pressures of the gas mixtures, blackbody target 
temperatures, instrument output voltages, and other housekeeping signals 
were recorded during each step of the calibration testing procedure. 

The recorded data were used in conjunction with a line-by-line radi- 
ative transfer program to compute the instrument radiometric calibration 
constants. For a given set of test conditions, this computer program calcu- 
lated the radiance values, L, aL, and AL', which corresponded to the test 
conditions for each of the measured signal output voltages of the three 
instrument channels, V, AV, and AV . 

The V channel calibration constants, R and S, for each instrument main- 
frame temperature, which is known to be the same as TB82, were determined by 
the least squares solutions to the calibration test data such that 

V » [R*(TBB2)][Lj3^g^j - Ljgj] * 6(TBB2) 

where r'^(TBB 2) is the responsivity (volts/watt cm-2 sr-^) and e(TBB2) is the 
offset voltage. The total inband radiances for the target blackbody, 
'‘target’ internal reference blackbody, Lggp were calculated 

using the line-by-line radiative transfer calculations. 
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The relationship between the output voltage signals for the av and AV' 
channels and the target radiometric signals AL and aL' is given by 


'"'target ’ a aV^^f(V.T6B2j 

'"'■target' t aVV,f(V.TBB2) 


Where R^'^(TBB2) 
and (TBB2) 


are the instrunent temperature dependent 
responsivities, 


and 


^Karget 

target 


are the calculated scene radiances for the AV and 
AV gas channels, respectively, and 


AV^g^[V,TBB2] 
and A VVef[''.TBB2] 


are the zero reference offset balance voltages for 
the AV and AV channels. 


The magnitude of the offset voltages is dependent both on the instrument 
temperature, TBB2, and on the magnitude of the target radiance as indicated 
by the V channel output signal, V. 

Preflight the MAPS instrument was very stable. The repeatability of 
the calculated constants was better than ± 2 percent for testing performed 
at approximately the same instrument temperature. The noise level on all 
three channels remained constant at 1.7 x 10^ watts • cm-^ * sr-^ for the V 
channel and 5.7 x 10-^ watts • cm-^ • sr-^ for the AV and AV channels, 

Postfl ight 

The large variation in coolant loop temperature encountered during the 
space flight necessitated an extensive set of postflight tests and 
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calibrations to characterize the instrument as a function of instrument 
temperature. The calibrations were carried out in a cryogenical ly pumped 
chamber capable of achieving a pressure of 1 x 10-'^ torr. The instrument 
was mounted on a heat transfer baseplate and suspended from the upper rail 
of the chamber. A secondary gas calibration reference cell and secondary 
standard apparatus replaced the MAPS flight camera. The gas cell was 
inserted and removed from the space between the optical port of the instru- 
ment and the secondary standard blackbody. 

The baseplate and secondary standard blackbody were connected to exter- 
nal thermal control units so that their temperatures could be varied inde- 
pendently. Thermistors similar to the internal thermistors of the MAPS 
experiment were used to measure the temperature of the secondary standard 
blackbody and the temperature of the 2.12-cm long secondary standard refer- 
ence gas cell. Test conditions are detailed in Table 1. 

The results from the entire set of test conditions were used to deter- 
mine the calibration constants for the instrument. The output of each chan- 
nel was an essentially linear function of the input signal. The respon- 
sivity, i.e., the slopes of the calibration curves, varied nearly linearly 
as a function of instrument temperature. However, the zero offsets of the 
AV and AV channels, in addition to being a function of the instrument tem- 
perature, were non-linearly related to the V channel signal. 

Infl ight 

As noted earlier, the MAPS instrument was exposed to very large coolant 
loop temperature variations during the STS-2 mission. These variations were 
so large and so rapid that it was not possible, in general, for the instru- 
ment to achieve a proper balanced condition during the 22-minute balance 
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Table 1. Test conditions for postflight instriinent characterization tests. 


Secondary Standard 

Instrument Temperature Secondary Standard Reference Cell 


(Kelvins) 

Target Temperature 
(Kelvins) 

Gas Mixtures 
(Percent CO by Volume) 

272.66 

273.16* 

2.5 

276.66 

280.16** 

3.8 

282.16 

293.16 

10.9 

289.16 

305.16 



312.16 



320.16** 



*Target temperature used in calibration of V channel only. 

**Target temperature equivalent to MAPS inflight calibration blackbody 
temperatures. 
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cycle. Only during the fourth rebalance attempt, which occurred about 13 
hours before the end of the mission, was a good balance achieved. 

The instrument temperature was almost the same during the fourth and 
fifth on-orbit balance/calibrate cycles, and it was nearly constant during 
the rebalance cycles; however, the instrument temperature varied between 
227 K and 285 K during the data acquisition period between these two bal- 
ance/calibrate cycles. Because the zero reference signal changes as the 
instrument temperature deviates from the balance temperature, it was neces- 
sary to determine the zero reference signal as a function of instrument 
temperature for the period between the fourth and fifth balance/calibrate 
cycles. This was accomplished by using data acquired while the instrument 
was viewing deep space, which is a known and stable source. 

STS-41G/0STA-3 Reflight 

In preparation for the reflight of the MAPS experiment as a component 
of the OSTA-3 payload, the instrument entered the T/V chamber again in Feb- 
ruary 1984. Appendix I contains the results of this calibration. 

The occurrence of the OSTA-1 temperature excursions made the MAPS ex- 
periment team aware of the need for additional inflight calibration data. 
Data takes with the MAPS sensor viewing deep space were added to the 41- 
G/OSTA-3 integrated payload command plan. These space look calibration 
points were taken at least twice daily during inertial maneuvers by the 
shuttle. Analysis of this inflight calibration is in progress. A better 
characterization of the zero reference signal as a function of instrument 
temperature will result. 

Two other significant changes to the MAPS experiment were made as a 
result of experience gained during the initial space flight of the MAPS 
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instriment. Electronic modifications were made to decrease the voltages of 
two signals, the vacuum channel signal and the housekeeping parameter R 
signal, to fall within the range of the A/D converter. Second, a tape 
recorder activate and inhibit command capability was implemented to enable 
selective data sampling. 

A transportable IBM-PC system was acquired for use by the OSTA-3 pay- 
load operations control team who worked at JSC. The potential exists for 
this computer system to provide some real-time data analysis. However, the 
system's primary functions during the OSTA-3 flight of the MAPS experiment 
were as monitor and recorder. Software modules were written to 

- display times and MAPS experiment status 

- enter comments into a flight log 

- print all or part of the flight log 

- simulate actual commands telemetered to the instrument 

- record onboard temperatures 

- display a directory of potential experiment malfunctions with 
corrective action(s) 

- alert experimenter to test conditions such as rebalance due or 
approaching end of recording capacity. 

Information acquired with this personal computing system has proven useful 

in the postflight OSTA-3 data reduction, which is now in progress. 

Data Reduction 

Major steps in the procedure whereby CO mixing ratios were inferred 
from the OSTA-1 MAPS experiment data are shown in Fig. 2. ODD researchers 
contributed in varying degrees to the successful completion of each branch 
of this data reduction method. In addition to work with the instriment 
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Procedure for the inferral of carbon monoxide mixing ratios. 







calibration as previously discussed, ODU personnel had primary responsi- 
bility for the lower branch of the flow chart. 

The MAPS experiment did not measure atmospheric temperature, so this 
data was obtained from the FNOC in the form of gridded meteorological north- 
ern and southern hemispheric analyses for 0 Z and 12 Z of the mission 
period. To incorporate atmospheric temperature into the regression method, 
the temperature profile was described as a single parameter, T^p, defined as 

P=0 . _ 

/ T(P) . SF{P)dP 

P=1000 ^ 
j _ n't* 

SF = p=o 

/ SF(P)dP 

P=1000 ^ 
mb 

where P = the pressure of the standard FNOC pressure levels 
T = the mean temperature between levels, and 
SF = the value of the signal function for each of the layers. 

The signal by the instrument received has the inherent property of being 
integrated over altitude. In a similar manner, the T^p is an average tem- 
perature, vertically weighted (by the signal function) over the depth of the 
atmosphere. Gridded T^p values were computed for each gas cell channel and 
for each analysis. Through spatial and temporal linear interpolation, the 
gridded T^p fields, in combination with the spacecraft ephemeris data, were 
used to produce values of T^p corresponding to the latitude, longitude, and 
time of each location where CO data were acquired. 

The character of the underlying surface and the measure of the solar 
zenith angle are two other parameters provided as part of the met file. The 
terrain data were selected from a Defense Mapping Agency terrain analysis of 
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1* X 1“ resolution. Solar zenith angles at each subsatellite point were 
calculated as a function of latitude, longitude, and time. 

DATA DISPLAY TECHNIQUES 

Computer generated graphics were used extensively as a method of sum- 
marizing, analyzing, and presenting data. Offsite, the hardware used for 
plotting was the Tektronix 4052 desktop graphics system. Onsite at LaRC, 
plot vector files were generated by a network of CDC Cyber mainframes, and 
images were postprocessed to drive the Varian, Versatec, and CalComp plot- 
ters. Both black on white and multicolored plots were produced. 

Mathematical and statistical software marketed by Tektronix was exe- 
cuted to generate hundreds of working plots like the preflight calibration 
curve shown in Fig. 3. These graphs provided quick looks at the linearity 
of data sets. The software allowed for data entry from keyboard, disk, or 
tape, and computed goodness of fit measures including R-square, residual 
error, and maximum absolute error. Both simple and polynomial regression 
curve fitting was available, and regression coefficients could be calculated 
and displayed. 

A general utility plotting program was coded in-house and executed 
extensively on the Tektronix system. The program displayed one or more sets 
of coordinate pairs on either linear or logarithmic axes. The user con- 
trolled scale factors, titles, synbols, and optional least squares fits to 
the data. Temperature profiles plotted with this utility are shown in Fig. 
4. The temperature values, taken from FNOC gridded analyses, were averaged 
within five latitudinal bands over the time period of the OSTA-1 flight. 

Such profiles were needed for the construction of the atmospheric models 
required in the data reduction process. 
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A graphics package originally developed at the National Center for 
Atmospheric Research (NCAR) has been implemented at the LaRC central com- 
puting complex. Among other capabilities, this package generates contour 
plots of equally spaced data and geographic outlines in nine different pro- 
jections. This software was used to produce Fig. 5. It shows 50 mb temper- 
ature fields taken from an NMC analysis and drawn on top of a polar stereo- 
graphic projection of the northern hemisphere. Contour plots of this type 
were used to compare meteorological analyses to actual measurements of atmo- 
spheric parameters. 

Figure 6 combines the NCAR graphics capabilities with the standard LaRC 
graphics library programs. It enables a rapid association of inferred CO 
values with the locations where the data were acquired. 

Another use of computer generated graphics as a data analysis tool is 
illustrated in Fig. 7. The cyclic variations with time of the V channel 
signal and of the internal reference blackbody (BBl) temperature are evi- 
dent. This and similar plots revealed the two second time delay in the V 
channel response to a change in the temperature of BBl. 

INFORMATION PROCESSING 

In parallel with the development of the MAPS experiment itself, the 
experiment team interfaced with other persons and organizations involved 
with the entire OSTA payload and the Space Transportation System. Through 
working group sessions, telecons, documentation, and memos, operational 
procedures were agreed upon which maximized the achievement of objectives by 
all experimenters and complied with shuttle operational and safety require- 
ments. ODD researchers participated in meetings such as the Payload Oper- 
ations Working Group (POWG) at JSC and the Payload Ground Operations Working 
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Fig. 5. 50 mb temperature fields for October 6, 1984, 12 Z. 
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Fig. 7. Cyclic variations with time of internal reference blackbody (BBl) temperature and V 
channel voltage. 







Group (PGOWG) at KSC. /Vnong the volumes of documentation which ODD person- 
nel helped draft, review, modify, and review again were the Test and Assem- 
bly Procedures (TAP), the Ground Interface Requirements Document (GIRD), and 
the Payload Integration Plan (PIP). 

ODD researchers were active in other avenues of scientific and techno- 
logical exchange. They gave oral presentations at events ranging from uni- 
versity seninars to international scientific conferences. They authored and 
co-authored technical papers for both internal use and public circulation. 
Appendix II is a listing of the journal articles, technical papers, and oral 
presentations of ODD personnel during the period covered by this report. 

The MAPS experiment has been guided since its inception by a Science 
Team drawn from academia and research organizations located across the 
nation and abroad. Not only have these scientists provided expertise during 
the design and data collecting phases of the experiment, but also they have 
interest in further specialized study of the data acquired. To facilitate 
their studies, members of the Science Team have been provided data ranging 
in format from copies of magnetic tapes containing meteorological analyses 
to multicolored plots of the global distribution of inferred CO. 

The OSTA-1 data set is currently available to the scientific community 
at large through the National Space Science Data Center (NSSDC) in 
Washington, D.C. In addition to the inferred CO data, relevant time, tra- 
jectory, and meteorological parameters, along with supporting documentation, 
were turned over to the NSSDC. Likewise, when the OSTA-3 data reduction has 
been completed, this data set will also be achieved at the NSSDA. 
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APPENDIX I 

□DURF TN NCCl -34-85- 15 

MEMORANDUK 

TO: Warren Hyoes, 

FROM: Joe Casas and Mary Saylor 

SUBJECT: MAPS OSTA-3 Calibration 

DATE: Auciust 5, 1985 


V CHANNEL 

1. Generate time history o-f the emitting sur-face o-f BB-1 

a. •£ period o-f BBl (seconds) is given bv a si mole 

regression eguation (See attachment A) 

b. ■£. peak to peak amplitude is 0,33 K and is independent 

o-f main-frame temperature TBB2 

c. V channel signal output lags the BBl thermistor readme 
by 2 seconds and is independent o-f TBB2 


2. N = (1/R) * (V - OFFSET) N 

DATA DATA BBl 

a. R = V channel responsivitv as a -function o-f TBB2 (K) 

(See attachment B) 

b. OFFSET = V channel o-f f set voltage as a -function o-f 

TBB2 (K) (See attachment C) 

c. NBBl = radiance of BBl as a -function of TBBl (K) 

(See attachment- D) 


and Qh#Q.Q®i.s 

1. Calculate Av and AV’ voltages using the 

zero ref zero ref 

bivariate fit coefficients (See attachments E and F) . 


a. AV (V.TBB2), where V is the V channel signal in 

zero ref 

volts and TBB2 is the mainframe temperature TBB2 (K) 


I-l 


b. 4V’ (V,TBB2>. where V is the V channel signal in 

zero re*f 

volts and TBB2 is the main-frame temoerature TBB2 (K) 

c. Prooram CKFITa4/UN*777475E 


a. dN = (1/R ) » ( 4V - ZV ) 

DATA ZV DATA zero re-f 

R * the resDonsivity o-f the ZV channel as a -function 
ZV 

o-f TBB2 (See attachment S) 


b. ZV-' 


DATA 


ZV’ 


(1/R ) » ( ZV’ - ZV ) 

ZV’ DATA zero re-f 

the resDonsivitv o-f the ZV’ channel as a 

-function o-f TBB2 (See attachment H) 


3. This approach should be used on the OSTA-3 -flight data during 
periods o-f temperature stable conditions -for TBB2, i.e., when 
the data are taken at a TBB2 temperature which is within 
1.5 K o-f the TBB2 temoerature during the balance cycle. 


4. Results o-f the ZV’ channel responsivity (spread in data) 
re-flect a need -for -further study as previously discussed. 
Several approaches could be taken. 
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EQUATION Y = A + B/K HAS MINIMUM MAXIMUM ABSOLUTE RESIDUAL 

SING SIMPLE REGRESSION USER DEFINABLE KEYS, SELECT THE FIT YOU WANT 
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